Astronomy & Astrophysics manuscript no. 5484tex 


© ESO 2008 


February 5, 2008 





Two-component magnetohydrodynamical outflows 
around young stellar objects 

Interplay between stellar magnetospheric winds and disc-driven jets 

! Z. Meliani 1 ' 2 3 , F. Casse 2 , and C. Sauty 1 

O ■ 

' Observatoire de Paris, L.U.Th., F-92190 Meudon, France 

u-p., e-mail: zakaria .meliani@obspm. fr, Christophe.Sauty@obspm.fr 

^ ■ 2 AstroParticule & Cosmologie (APC)* - Universite Paris 7, 11 place Marcelin Berthelot, F-75231 Paris Cedex 05, France 

' e-mail: fcasse@apc.univ-paris7.fr 

3 Max Planck Institute for Astrophysics, Box 1317, D-85741 Garching, Germany 



C^*) , Received ... / accepted ... 

ABSTRACT 



(N 
> 

in 
o 



Context. We present the first-ever simulations of non-ideal magnetohydrodynamical (MHD) stellar magnetospheric winds coupled with disc- 
driven jets where the resistive and viscous accretion disc is self-consistently described. 
qq , Aims. These innovative MHD simulations are devoted to the study of the interplay between a stellar wind (having different ejection mass rates) 

" and an MHD disc-driven jet embedding the stellar wind. 
\£) Methods. The transmagnetosonic, collimated MHD outflows are investigated numerically using the VAC code. We first investigate the various 
f*"*) . angular momentum transports occurring in the magneto-viscous accretion disc. We then analyze the modifications induced by the interaction 
between the two components of the outflow. 

Results. Our simulations show that the inner outflow is accelerated from the central object's hot corona thanks to both the thermal pressure 
and the Lorentz force. In our framework, the thermal acceleration is sustained by the heating produced by the dissipated magnetic energy due 
O ' to the turbulence. Conversely, the outflow launched from the resistive accretion disc is mainly accelerated by the magneto-centrifugal force. 
4^3 Conclusions. The simulations show that the MHD disc-driven outflow extracts angular momentum more efficiently than do viscous effects in 
ffi ■ near-equipartition, thin-magnetized discs where turbulence is fully developed. We also show that, when a dense inner stellar wind occurs, the 
. . resulting disc-driven jet has a different structure, namely a magnetic structure where poloidal magnetic field lines are more inclined because of 
the pressure caused by the stellar wind. This modification leads to both an enhanced mass-ejection rate in the disc-driven jet and a larger radial 
extension that is in better agreement with the observations, besides being more consistent. 
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1. Introduction 

Accreting stellar objects are often associated with collimated 
jets or winds from accretion discs. Most of those objects also 
show evidence of winds originating in a corona surrounding 
the central object. These accretion-ejection phenomena are ob- 
served in different astrophysical sources ranging from young 
stellar objects (YSOs), X-ray binaries, planet ary neb ulae to ac- 
tive galactic nuclei (AGNs) (see, e.g. dLiviolll997h and refer- 
ences therein). The outflow provides an efficient way of ex- 
tracting angular momentum and converting gravitational en- 
ergy from the accretion disc or from the central object into out- 
flow kinetic energy. 
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Observations show that most of the jets are launched very 

close to the central engi ne. In the case of YS Os, there is direct 

J II "i 

observational evidence (Burrows et al 1996) as well as in the 

case of some microquasars ( Fender et allll997tlMirabeill2003l) . 
For instance, it has been suggested that, in microquasars, the 
fastest components o f the o utflow are launched in the vicin- 
ity of the black hole jMeien lt2003). Another piece of evidence 
that the outflow may originate in a region relatively close to the 
central object is that the observed asymptotic velocity of the 
jet is close to the escaping speed from the central engine. Thus 
there is a direct relation between the asymptotic speed and the 
depth of the g ravitational potential ( Mirabe ll 19991 iLiviol 1 999: 
Pringle 1993). Moreover, the high-resolution images of H a and 
TO/I dBacciotti et all2002l) reveal a continuous transverse vari- 
ation of the jet's velocity, where the fastest and densest compo- 
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nents are closer to the central axis. 

The high velocity of the observed jet in YSOs suggests that 
they originate in a region th at is no larger than one a stronom- 
ical unity (AU) in extent jKwan & Tademarulll988l) and be- 
tween 0.3 to 4.0 AU from the star in the case of the LVC of 
DG Tau dAnderson et allE003l) . This theoretical prediction is 
supported, in the case of a disc wind, by observ ations of the 
rotati on of several jets associated with TTauris dCoffev et all 
2004- Moreo ver, in the case of classical TTauris (CTT S), UV 
observations ( Beristain etail200lt lOunree et all 2005) reveal 
the presence of a warm wind whose temperature is at least of 
3 x 10 5 K. It appears that the source of this wind is restricted to 
the star itself, where X-ray ob servations support the presenc e 
of a hot corona in CTT stars (Feigels on & Montmerld fT999). 
These observations also suggest the existence in CTTS of stel- 
lar winds comparable to the solar wind. These winds may be 
both thermally and magneto-centrifugally accelerated. 

Since the discovery of the existence of winds and jets in 
astrophysics, enormous progress has been made regarding the 
understanding of these phenomena. At the same time, we still 
do not know precisely how the wind from the central corona of 
the star or the compact object interacts with the disc outflow 
and the respective roles and differences between these two 
types of flows. 

Several works have analytically and numerically stud- 
ied the formatio n of outflows la unched from the ac- 
cretion disc i jBJandforfjfePavne] 
IContopoulos & Lovelace! 



1994 



1995 



199 



lOuved et all Il997t IVlahakis & Tsingano; 



Casse & Ferreiral l2000at ICasse & Keppensl I2C 
Anderson et al 2005; Pudrit z et all 1200(1) . Other works have 




focused on the outflows from the hot corona of the c entral 
objects llSakurailll985USautv & Tsinganoslll994UFendtll2003t 
iMatt & Balic kl2004 . 

In models dealing with outflows launched from accretion discs, 
the magnetic field plays a key role in the accretion, the accel- 
eration, and the collimation of the associated v ertical wind, 
which is also supported by recent observations iDonati et all 
2005). 

The detection of rotation signatures in TTauri jets gives 
extra strong support to the magneto-centrifugal launching 
from the accretion disc. However in stellar ouflows, the wind 
must be thermally accelerated because of the strong heat- 
ing of viscous and non-ideal magnetohydrodynamical (MHD) 
mechanisms. This acceleration increases and approaches the 
magneto-centrifugal acceleration at least near the polar axis. 
Some models have already investigated diffusive disc-driven 
jet launching. In some simulations the accretion disc was con- 
sidered as a fixed, time-independent boundary condition, while 
a con stant magnetic resistivity prevails through the e ntire out- 
flow (fFendt & Cemeliicl2002l) . lKuwabara etal d2005t) have in- 
cluded an accretion disc in their resistive simulations but only 
considered an uniform magnetic resistivity everywhere in their 
computational domain. Because of the very different physi- 
cal conditions prevailing in the disk, the jet, and the external 
medium, it seems very unlikely that this kind of resistivity is 
physically relevant. Models involving two component bipolar 
outflows have been proposed in the case of AGN as such as 



in Sol et al. (1989) and Renaud & Henri (1998), where an 
electron-positron central wind component is surrounded by an 
external, ideal MHD disc-driven jet. Another two-component 
outflo w model has been proposed by ( Bog ovalov & T singanos 
2005) where a relativistic pulsar wind is surrounded and colli- 
mated by an ideal MHD disc-driven wind. In the case of YSOs, 
two-comp onent models were considered, such as in X -wind 
outflows dSautv & Tsinganos! Il99l iFerreira et alll2000l) . The 
inner component extracts its energies from the corona around 
the central region (the central object and the inner part of the ac- 
cretion disc where an advection-dominated accretion flow may 
exist), while the second component is launched from the thin 
accretion disc. 

The aim of this paper is to investigate the formation of a two- 
component outflow around YSOs, one coming from the thin 
accretion disc and the other one injected from the hot corona of 
the central star. This work is developed on the base of the disc 
wind simulations of Casse & Keppens (2002, 2004) (CK02, 
CK04). Its goal is to study the influence of the stellar wind 
on the structure and the dynamics of the jet around YSOs. 
Furthermore, we investigate the consequences of the energy 
dissipation in the outflow close to the polar axis. Before that, 
we present simulations of the outflow launched from a resis- 
tive and viscous accretion disc where magnetic Prandtl number 
(ratio of the anomalous viscosity to the anomalous resistivity) 
equals unity. This is the first time that viscosity and resistivity 
are considered together in the disc and included in an MHD 
simulation involving both the accretion disc and the related 
jet. In previous works, the viscous accretion disc was exam- 
ined without taking the disc wind into account or with an im- 
posed internal structure rtVon Rekowski et all2006l) . which does 
not enable the authors to study the complete angular momen- 
tum transfer. Thus, the first part investigates the relative role of 
angular-momentum transport by viscosity and by the outflow- 
ing plasma and its influence on the formation of the outflow 
from the disc. Then in a second part, we present the results 
of our simulations of ideal MHD outflows launched from re- 
sistive, viscous, accretion discs surrounding the turbulent wind 
accelerated from the hot corona of the central star. 



2. Ideal MHD outflows arising from resistive, 
viscous, thin accretion discs 

In this section we present the simulations of axisymmetric 
MHD outflows generated from thin accretion discs where for 
the first time viscosity is taken into account, together with 
resistivity. We recall that turbulence is believed to generate 
both anomalous resistivity and viscosity, such that the turbu- 
lent magnetic Prandtl number, which is the ratio of the viscous 
to the resistive transport coefficients, should be of order unity 
within flows where turbulence is f ully developed ( Pouau et et al 
1 1974 Kitchatinov & Pipin 1994). This is supposed to be the 
case of accretion discs and stellar winds. The presence of two 
braking torques inside the disc may achieve different disc-jet 
configurations, since the angular momentum transport is mod- 
ified by the presence of viscosity. 
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2.1. MHD equations 

To get the evolution of such a disc, we solve the system of 
time-dependent resistive and viscous MHD equations, namely, 
the conservation of mass, 



and momentum, 

dpv (B 2 

+ V (vpv - BB) + V — + P + pVO G = 
at \ 2 ' 



(1) 



(2) 



We also consider the energy conservation governing the tem- 
poral evolution of the total energy density <?, 



B 2 pv 2 P 
2 2 y - 1 



(3) 



<3e 



+ V' 



BB v 



i] m J 2 - B ■ (V x r, m J) - V(v?7 v ft) 

where p is the plasma density, v the velocity, P the thermal pres- 
sure, B magnetic field, y the specific heat ratio (Cp/Cy = 5/3), 
and J = V x B is the current density. In this set of MHD equa- 
tions, the thermal pressure is derived from all conserved phys- 
ical quantities. To close the set of equation, we have included 
a perfect gas equation of state linking the thermal pressure and 
plasma density to the plasma temperature T, namely T - P/p. 
The gravitational potential is given by 




Fig. 1. Density contours in the poloidal plane of an accretion- 
ejection structure where a viscous and resistive MHD disc is 
launching a collimated jet. Magnetic field lines are drawn in 
black solid lines, while the fast magnetosonic surface corre- 
sponds to the white solid line (Alfven surface is the black dot- 
ted line). The size of the sink region is = 0.1AU and the 
stellar mass is 1M . 
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Note that both resistivity (t],„) and viscosity (77,.) are taken into 

(R 2 + Rlf 

account in the MHD set of equations. The viscous stress ten- £ 

sor is given by rj v fl = -rj v ((Vv + Vv T ) + |/(V • v)). The last V g (R,Z) = V (l - e 2 ) 2 _ exp 

equation provides the temporal evolution of the magnetic field, e ^2 + ^2^ 4 

namely the induction equation 



R 
V 7 - 



H 2 I ~Z 
2Z 2 ^ 



H 2 



(7) 



(8) 



dB 

— +V(vB-Bv) = -Vx(r ]m J) . 
at 



(5) 



The local heating and torque in the accretion disc are gener- 
ated by the magnetic resistivity and hydrodynamical viscos- 
ity occurring in the disc. We adopt a magnetic Prandtl number 
Pr — — — 1 in our simulations as it is reasonable on physical 
grounds (see references above). 



2.2. Initial and boundary conditions 

The initial density profile and velocities are set as follows, 



where H = eR is the disc height, which is proportional to the 
radius R, via the disc aspect ratio e ~ C s /Vk linking the disc 
sound speed C s to the Keplerian velocity Vk- We deliberately 
choose £ smaller than unity to get a thin disc where thermal 
pressure gradi ent is smaller than both centrifugal and gravita- 
tional forces ( Wardle & Konial 1993). The parameter m s = 0.1 
ensures that the initial poloidal flow remains subsonic. 
The initial magnetic field configuration is taken as in CK04, 



F{R,Z) = ^,B 



Rl /4 R 2 



1 



B R (R,Z) = - 



1 dF 
Y&Z ' 



p(R, Z) = p max 



5 x 10 



-6 



R- 



(R 2 + RlY 



B Z {R,Z) 
B g (R,Z) 



V06z) 



1 dF 
RdR ' (1 +R 2 ) ' 
0. 



(9) 



5xlO~ 6 ,<M - 



(r-i)z 2 



H 2 



(6) 



where /3 P = B 2 /P is set to 0.6 to ensure that the magnetic 
pressure remain on the order of the thermal pressure, a nec- 
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essary condit ion for a dis c launching large-scale stable jets 
( Ferreira & Pelletier|[l995l) . The parameter £ controls the ini- 
tial bending of the magnetic surface. In the following simula- 
tion we set f = 0.04. The reader may shift from dimensionless 
quantities used in our simulations to physical quantities by set- 
ting the mass of the star M», the accretion rate M a , as well as 
the size of the disc inner radius /?;, 
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The expression of the anomalous resistivity in the accretion 
disc is as in CK04, 



rjv_ 
Pr 



a m V A \z=oHexp[-2 



H 2 



(11) 



The anomalous resistivity r\ m has the same origin as the tur- 
bulent viscosity since time correlations for hydrodynamic and 
magnetic turbulence are the same. Thus we likewise intro- 
duce an anomalous viscosity r\ v equal to r\ m . As in CK04, we 
have replaced in the a-prescription the sound speed by the 
Alfven speed because the anomalous resistivity is related to 
the small-scale turbulent magnetic field. Since the disc remains 
near equipartition between thermal pressure and magnetic pres- 
sure, this does not affect the value of the transport coefficients. 
Through the dependence on the Alfven velocity, this becomes a 
profile varying in time and space that essentially vanishes out- 
side the disc. We take a m = 0.1 smaller than one to ensure that 
the Ohmic dissipation rate at the mid-plane of the accret ion disc 
does n ot exceed the rate of gravitational energy release ( Konigl 
1 19951) . However, for a lower value of a m < 0.1, the energy 
released by accretion is insufficient to produce a strong colli- 
mated wind crossing all critical surfaces. In the case of weak 
resistivity, the resulting o utflow re mains weak, and the opening 
angle of the jet is small dFerreiralll997l). This result has bee n 
confirmed by numerical calculations dCasse & R eopens 2002). 
The boundary conditions here are similar to CK04. We de- 
signed an absorbing sink around the origin in order to avoid the 
gravitational singularity. In the first quadrant of the simulation, 
the sink region is a square of one unit length both in the R and 
Z directions where matter can only enter the zone {Vr, Vz < 1) 
in order to avoid any numerical artifact. We consider the axis 
and disc mid-plane as a combination of symmetric and anti- 
symmetric boundaries. The top and right boundaries are set as 
free boundaries (with nil gradients) except for the outer radius 
of the disc where we impose a fixed poloidal mass accretion 
rate. 

The numerical simulations presented in this pap er were per - 
formed using the Versatile Advection Code VAC ( lT6thll9 96). 
see http://www.phys.uu.nl/~toth We solve the full set 
of resistive and viscous MHD equations under the assumption 



Fig. 2. Temporal evolution of several angular momentum fluxes 
occurring inside the accretion-ejection structure displayed in 
Figffl The various fluxes are normalized to the amount of an- 
gular momentum removed from the disc Lmec- The dominant 
way to remove disc angular momentum is provided by the mag- 
netic torque leading to the creation of a jet. 




Fig. 3. Same as in Fig^but with a ideal stellar wind emit- 
ted from the inner region whose ejection mass rate is M = 
lO~ 9 M Q /yr. The disc-driven jet conserves a very similar dy- 
namical structure to the case where no stellar wind is emitted. 
The size of the sink region is Ri = 0.1 AU and the stellar mass 
is 1M . 



of cylindrical symmetry. We time-advance the initial conditions 

using the conservative, second-order accurate total variation di- 

i II 
minishing Lax-Friedrichs scheme (Toth & Odstrcil 1996) with 

minmod limiting applied to the primitive variables. We apply a 
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projectio n scheme prior to each time step to enforce V • B — 
jBrackbill & Barnesll980l) . 

2.3. Angular momentum transport in resistive, viscous 
thin accretion discs launching MHD jets 

In order to study the angular momentum transport governing 
such an accretion disc, we direct our first simulation here to 
studying the sole resistive and viscous accretion disc threaded 
by a large-scale magnetic field. We do not set any additional 
outflow coming from the central star in this simulation. 
Writing down the conservation of angular momentum in an ax- 
isymmetric framework, we see the two mechanisms responsi- 
ble for the angular momentum transport and removal, namely 
the magnetic torque and the viscous torque 

dpVg 



dt 



+ V ■ (pRVgV p - RBgB p + Rpj] v n p9 ) = 



(12) 



where the subscript "p" stands for the poloidal component of 
the labeled vectors. These two toroidal forces allowfor angu- 
lar momentum transport in two different directions. Indeed, 
the magnetic torque provides angular momentum along the 
poloidal magnetic field, namely in the jet direction. Since thin 
accretion-disc plasma velocity is dominated by the Keplerian 
rotation, the viscous stres s tensor can be approximated as 
JShakura & Sunvaevll973l> 



Rrjvpllpe =s r] v Rp^^-e R = ~pr] v RQ. K e R 



(13) 



where is the Keplerian angular velocity. The viscosity will 
then initiate a radial angular momentum transport. The range 
of turbulence configuration in the accretion disc is quite end- 
less (as the range of Prandtl number), so we restrict ourselves 
to the configu ration predicted in the case of fully develop ed tur- 
bulence (Poua uet et a! 19761 iKkchatinov & Pipinll994l) . 
The result of the simulation is displayed on Fig^ logarith- 
mic density; the poloidal magnetic field lines. The obtained 
accretion-ejection structure shows a super-fastmagnetosonic 
collimated jet. In order to study the difference with the 
accretion-ejection flows obtained by CK04, where the mag- 
netic Prandtl number was set to zero, we measured the var- 
ious angular momentum fluxes crossing the internal and ex- 
ternal radii, as well as through the disc surface. To do so, we 
define global variables characterizing the angular momentum 
extracted from the accretion disc, namely, 



'-'LIB - 

where 

Lmec 



^MEC 



+ L 



MHD 



+ L 



vis 



(14) 



= - ff dS lP vRV 9 - ff aS E -pvRV e (15) 
JJs, JJs E 

is the variation of the angular momentum advected by the in- 
flow between the external and the internal parts of disc, 



Lmhd 



ff dSiB p RBg - ff dS E ■ BpRBg 
JJs, JJs E 



(16) 



is the variation of the MHD Poynting flux between the internal 
and external radii. This magnetic contribution to the radial an- 
gular momentum accounts for the twisting of the magnetic field 



occurring inside the disc, which is similar to storing angular 
momentum and mechanical energy of the plasma in the mag- 
netic field (generating toroidal magnetic field Bg). The amount 
of disc angular momentum removed by viscosity is given by 



^vis 



ffs. 



dSi 



■ eRpj] v RYl R g - dS E ■ e R pj] v RYl R g. (17) 



We denoted by dS E = 2nR E dZeR the outer and by dSj = 
InRidZeR the inner vertical cut through the accretion disc, with 
-H < Z < H. 

When the outflow is arising from the accretion disc, we can 
evaluate the angular momentum transported vertically into the 
jet by considering the various fluxes through the disc surface, 



+ Ly 



where 

Lmecj 



dS surf ■ pV p rVg 



(18) 



(19) 



is the angular momentum advected by the vertical mass flowing 
into the jet, 



dS sur f ■ BpRBg 



(20) 



is the angular momentum extracted by the magnetic torque 
in the accretion disc and converted into MHD Poynting flux 
through the disc surface. Finally, 



^VISJ 



ffs s „ lf 



dS surf ■ e z pri v RIl z g 



(21) 



is the angular momentum extracted by the viscous torque at the 
disc surface where r} v YV z g = -t] v R'-^. However, the effect of this 
mechanism is zero, as viscosity vanishes outside the accretion 
disc (Eq. lilt . We note the disc surface as dS sur f = InRdRez 
with Ri < R < Re- In the case where the whole structure 
reaches a stationary state, the angular momentum conservation 
can be translated into a global relation Lub - Ljet- The angu- 
lar momentum fluxes related to our simulation are displayed in 
Fig|5] where we represent L MHD , L VIS and L JE t normalized to 
the flux of angular momentum removed from the disc L ME q. 
It clearly shows that inside a resistive, viscous, thin accretion 
disc with Prandtl number equal to unity, the viscosity is unable 
to remove the disc angular momentum efficiently (as already 
showed by Pudritz & Norman 1986), since only one percent is 
carried away by the viscous torque. Conversely the presence 
of the jet has an strong impact on the angular momentum 
balance because it enables the magnetic torque to achieve a 
very efficient angular momentum transport from the disc into 
the jet (more than 90%). Among the three components of 
Ljet, L vis j is nul because viscosity is vanishing at the disc 
surface and L MEC J is small compared to the MHD Poynting jet 
flux Lmhdj (mainly because mass is sub-slow-magnetosonic 
at the disc surface). This magnetic energy reservoir created at 
the base of the jet is used in the jet to accelerate matter such 
that the jet becomes super-fast-magnetosonic. It is noteworthy 
that the structure coming from our simulation reaches a 
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Fig. 4. Time evolution of a two-component jet launched from a thin accretion disc threaded by a bipolar magnetic field. The 
outflow is composed of a disc-driven jet embedding a non-ideal stellar wind emitted from a YSO located at the center of the 
simulation in the sink region. The density contours are represented by greyscales while poloidal magnetic field lines are displayed 
using solid lines. The various snapshots represent the same system but at different stages (after five, ten, twenty, and thirty 
inner disc rotations). The simulation was performed with a stellar mass-loss rate of M = lO~ 9 M G /yr. In the early stages of our 
simulation the stellar wind is the only outflow of the system while as the simulation goes on, a disc-driven jet appears around the 
stellar outflow and collimates it. The size of the sink region is R/ = 0.1 AU and the stellar mass is 1M Q . 



quasi-stationary state where Ljet — Lub- 

In order to explain why the magnetic torque prevails in the 
accretion disc, we can use some ordering to estimate the rela- 
tive amplitude of the two torques. The magnetic torque expres- 
sion is 



/ f v n\ „dB e B R dRB e 

(J X «) ■ e e = B z —— h — 

dZ R dR 



(22) 



that we can simplify knowing that Br <k Bz inside the disc and 
that the toroidal magnetic field that is nil at the disc mid-plane 
approaches —Bz on the disc surface, 



B 2 

|(/xB).e^-|. 

ti 



(23) 



The viscous torque, as previously written, can be expressed as 



|V • (flj/.Upe ■ e e )\ - 



dpr] v RCl K 



dR 



R 



(24) 



since the radial variatio n of the various qu antities is expected 
to be like power laws iShakura & Sunvae\ll973h . The ratio of 
the two torques can be written as 



( J x B) ■ e e 



(V ■ Rr] v n p9 ) ■ e 9 



Va 
Cs 



1 

a v e 



(25) 



where the disc aspect ratio e is much lower than unity in a 
thin accretion disc, as is the viscosity parameter a v . The ac- 
cretion disc has to be close to equipartition between mag- 
netic pressure and therma l pressure in order to launch a jet 
( Ferreira & Pelletien l 19951) . so it is easy to see that this ratio 
is much higher than unity in all magneto-viscous, thin-disc 
launching jets. In the context of our simulation, we have set 
e = a v = 0.1 leading to a ratio on the order of 10 2 , which 
is compatible with the ratio of L jet to Lyis in Fig0 In con- 
clusion, we see that angular momentum extraction from a thin 
or even a slim magnetized-disc (e <K 1) is likely to occur 
in the disc-driven jet rather than in the disc itself, for a disc 
close to equipartition, i.e. with a plasma beta close to one. This 
results is consistent with previous analytical models of non- 
resistive disc winds where the accreti on-related wind removes 
the excess of the angular momentum ( Pudritz & Norman 1986; 
IPelletier & Pudritzll992tlLubow et al. Il994 . 
However, in order to have a more consistent simulation of the 
accretion-ejection structures, we should take the interaction 
with the inner stellar coronal wind into account. Incontrast to 
previous simulations, we include the acceleration of the stellar 
wind, which is likely to start with a subsonic motion from the 
base of the corona and then accelerates, as well as the full de- 
scription of the accretion-disc launching jets. The stellar wind 
acceleration close to the axis cannot be exclusively magnetic, 
since magneto-centrifugal effects vanish near the axis. Besides, 
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Fig. 5. Left: plot of the initial temperature isocontours of the accretion disc. Right: plot of the temperature isocontours of the 
simulation displayed in Fig[3] Temperature isocontours are not displayed in the external medium (outside of both the jet, stellar 
wind, and disc) as it is considered a near vacuum medium with very low temperature. The size of the sink region is Rj = 0.1 AU 
and the stellar mass is 1M . 



the high coronal temperature is likely to induce a more efficient 
turbulent heating. We thus intend to use the turbulent wind vis- 
cosity and resistivity as the primary sources of acceleration for 
the inner stellar outflow. Turbulence may be induced in the 
stellar magnetospheric wind by its interaction with the disc- 
driven jet. The differences between both their dynamics and 
thermodynamics probably induce instabilities. The turbulence 
may also have a stellar origin and/or a possible connection to 
the accretion occurring near the surface of the star. In fact, the 
inner accretion surface, as well as the star surface, are time- 
dependent and inhomogeneous, leading to outwardly propagat- 
ing Alfven waves in the stellar wind and inducing turbulence. 
This origin of turbulence is based on an analogy with mod- 
els and observations of the so lar wind where the solar o rigin 
of tur bulence is investigated (Leam on et all Il998l ISmith et a] 
l2000l) as the convection below the photosphe re (Cranmer & 
Ball egooijen 2005; ( See also a review papers dGoldstein et all 
ll995UCranmerl2004 l. 

3. Two-component MHD outflows from a resistive, 
viscous, accretion disc and a star corona 

The aim of this section is to model the interplay between the 
two components of a YSO outflow, namely a jet launched from 
a magnetized resistive, viscous accretion disc and the second 



one, an non-ideal MHD spherical wind ejected from the pro- 
tostar hot corona. In this section, we describe the method that 
we developed for the implementation of the stellar wind in the 
model. The results and the differences, with the non-ideal MHD 
stellar wind model, is discussed in the following section. We 
use the same initial conditions as in the previous section, ex- 
cept for is a change in the boundary condition located at the 
top of the sink region. 

We now replace the accretion inflow by an outward mass 
flux whose amplitude is r times the solar mass-loss rate and 
which is spherically ejected with a speed that is a fraction 5 of 
the fast-magnetosonic speed in the corona V"; 

M = Tl0- l4 M o /yr (26) 

v* = <5,V;| (27) 

v" z p = 6 v Vf. (28) 

In our computational domain, the sink region around the origin 
is a square of one unit length in both the R and Z directions. 
This length corresponds to the internal radius Rj of the accre- 
tion disc but does not correspond to the physical radial disc 
edge that is located inside the sink. In the following plots, we 
set Rj to 0. IAU so that the upper limit of the sink region and our 
internal disk radius are at 20 star radii from the star surface in 



8 



Meliani, Casse & Sauty: Two-component MHD outflows around YSO 




10 20 30 10 20 30 10 20 30 

Inner radius rotation Inner radius rotation Inner radius rotation 



Fig. 7. a) Plot of the ejection mass rate from the accretion disc, b) the inner accretion rate, c) the ratio of the stellar mass-loss rate 
to the disc-wind mass-loss rate as a function of time. These plots are related to the simulation performed with a stellar mass loss 
of M = KT 9 M B /yr. 
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Fig. 6. Plot of the vertical variation at R = 1 AU from the axis of 
the temperature in the simulation displayed in Fig|3] The size 
of the sink region is Ri = 0.1 AU and the stellar mass is IM@. 



the case of a one solar-mass central object. The magnetic field 
near the star has a near spherical expansion that is becoming a 
near vertical structure (Eq|9j. Our magnetic field prescription 
is coherent with a magnetic field at the surface of the star on 
the order of B ~ 2kG in agreement with o bservational values 
Jjohns-Krull et afcOOlllGuentther et al 19991) . 

In our simulations, we model neither the very inner part of 
the accretion disc nor its interaction with the magnetosphere of 
the star located at (3R Q - 6 R Q ). 

Nevertheless, we include the effect of the star rotation in our 
simulations by imposing a solid rotation profile (Vg/R = est) on 
the outflow at the top of the sink region. We set the angular ve- 
locity of the outflow at this boundary to the Keplerian angular 
velocity at the inner radius Ri. The rotation period associated is 
then 

/ r \ 3/2 /m r 1/2 
p - = n - 57da HoIXu) fe) • (29) 



which corresponds to young star rotation periods such as GM 
Tau JGullburing et all998l) . 

We also consider two different cases of stellar winds. The 
first one is consistent with a heavy hot wind whose mass loss is 
on the order of M = lO~ 7 M Q /yr (r = 10 7 ). This mass-ejection 
rate is in the range of typical mass losses for young B-star 
and O-star type. This kind of YSO is characterized by strong 
outflows and dynamical timescales around I0 4 years, and the 
stellar wind is believed to be the main contributor to the out- 
flow. The computational domain related to this simulation is 
[R, Z] = [0, 80] x [0, 120] with a resolution of 304 x 404 cells. 
The other simulation stands for systems where the stellar wind 
is a light and hot one, namely with a mass loss on the order of 
M - lO~ 9 M /yr (t = 10 5 ) (FigGJ. In these cases the stellar 
wind is lighter than the jet launched from the accretion disc. 
The computational domain associated with the second simula- 
tion is [R,Z] = [0,40] x [0, 80] with a resolution of 134 x 204 
cells. It is noteworthy that, since we are expecting a widening 
of the jet due to the strong stellar mass loss in the first simula- 
tion, we have designed a larger computational domain in order 
to capture all the features of the resulting outflow. We set the 
velocity parameter 8 V to 0.01 in both simulations, which is con- 
sistent with an initial sub fast-magnetosonic and sub-Alfvenic 
ejection from the corona. We assume that part of the outflow 
acceleration has already taken place between the star surface 
(hidden in the sink) and the top boundary of the sink region 
although the flow is sometimes even sub-slow-magnetosonic, 
depending on the magnetic configuration. 

3.1. Non-ideal effects in stellar winds 

In most stellar wind models, the wind material is often subject 
to coronal heating, contributing to the global acceleration of the 
flow. In our simulations, we assume that the coronal heating is 
a fraction 6 e of the energy released in the accretion disc at the 
boundary of the sink region, which is transformed into thermal 
energy in the stellar corona close to the polar axis. This see- 
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nario was proposed by jMatt & Pudrit J2005b and is supported 
by the current observations of hot stellar outflows ( Dupre e~et all 
2005). The thermal energy imposed at the lower boundary of 
the corona is, at each step of the simulation, the sum of the 
thermal energy s up+l of the above stellar jet, i.e. the thermal 
energy of the first cells above the sink border, plus a fraction 
6 e of the thermal energy at the disc inner radius e'". Thus the 
thermal energy at the upper boundary of the sink is 



s = s up+i + s in 6 F . 



(30) 



The 5 e parameter range is limited from below by the initial 
thermal acceleration at the surface of the corona, which should 
balance the gravitational force, and from above by the require- 
ment of avoiding too high a temperature in the corona. In our 
simulation we take a small efficient heating corona S e = 10 1 . 
We deliberately use a very low value for this parameter in or- 
der to ensure that the main heating source of the stellar wind 
lies in the Ohmic heating. We intend to study its effects com- 
pared to the prescription of a larger amount of thermal energy at 
the base of the stellar flow. The low value of 5 e represents the 
amount of energy released by Ohmic heating below the sur- 
face of the sink. This heating is essential for letting the flow 
escape from the gravity, since the flow has already undergone 
an initial acceleration from the surface of the star to the top 
of the sink region. Moreover, the wind undergoes a mechani- 
cal heating where the accreted flow at the top of the accretion 
disc compresses the inner wind and may sustain turbulence in 
the wind. In Sect. 3.2.2, we will discuss the ideal stellar winds 
emitted from the sink where a large amount of thermal energy 
is deposited at the base of the flow. 

The interaction between the different components of the out- 
flow may be responsible for energy dissipation inside the 
plasma, which is the outcome of non-ideal MHD mechanisms 
occurring in the wind. In this paragraph, we show how these 
non-ideal MHD effects are taken into account when prescrib- 
ing a turbulent magnetic resistivity taking place in the wind 
region, in addition to the disc resistivity 

T] m = a m Va\ z =o H exp ( -2 ^ J 



+ a w V A H w exp 



R 



(31) 



The first term accounts for the anomalous resistivity occurring 
in the accretion disc. It vanishes outside the disc (Z > H). The 
second term corresponds to the description of an anomalous 
resistivity occurring in the outflow close to its polar axis. This 
term vanishes outside the stellar wind (R > H w ) where H w is the 
distance from the polar axis where the Alfven speed encounters 
a minimum. Hence, the dissipation effects are only located in 
the stellar wind component and not in the disc wind, which is 
supposed to be less turbulent. For the resistivity in the stellar 
wind, we take a w = 10~ 2 , a lower value than in the disc itself. 

3.2. Two-component MHD outflows from YSO 

We first focus on a simulation where the stellar mass loss is set 
to 10 _9 M o /yr. The outcome of our simulation can be seen in 



Fig.0]where we display four different snapshots of the poloidal 
cross-sections of the structure at different times its evolution, 
i.e. at 8, 16, 24, and 32 rotation periods of the inner disc ra- 
dius. In these snapshots we display the density contours and 
the poloidal magnetic field lines. The initial accretion disc con- 
figuration is close to a hydrostatic equilibrium where the cen- 
trifugal force and the total pressure gradient balance the grav- 
ity. In the central region, the matter is continuously emitted at 
the surface of the sink region (designed to be close to the star 
surface) with sub fast-magnetosonic speed and with a solid ro- 
tation velocity profile. Initially, a conical hot outflow (stellar 
wind) propagates above the inner part of the disc. Its inertia 
compresses the magnetic field anchored to the accretion disc. 
As a result, the bending of the magnetic surfaces increases, 
leading to a magnetic pinching of the disc. This pinching delays 
the jet launching as the disc has to find a new vertical equilib- 
rium. Thus the disc takes a few more inner disc rotations be- 
fore launching its jet compared to CK04. Once the jet has been 
launched, the structure reaches a quasi steady-state where the 
outflow becomes parallel to the poloidal magnetic field that is 
parallel to the vertical direction. 

The obtained solution is fully consistent with an accretion-disc 
launching plasma with a sub-slowmagnetosonic velocity. The 
solution crosses the three critical surfaces, namely the slow- 
magnetosonic, the Alfven, and the fast-magnetosonic surfaces. 
The other component of the outflow, namely the stellar wind, 
is injected with sub-fastmagnetosonic velocity and crosses the 
Alfven and fastmagnetosonic surfaces. The two components 
of the outflow become super-fastmagnetosonic before reaching 
the upper boundary limit of the computational domain. Figure^] 
also shows that the outflow has achieved quite a good colli- 
mation within our computational domain. We can distinguish 
between the two components using the isocontours of tempera- 
ture, which are displayed as grey-scales in Fig.|5] In this figure, 
we can clearly see a hot outflow coming from the central ob- 
ject embedded in the cooler jet arising from the accretion disc. 
In Fig.|5]we also show that the thermal energy released by the 
Ohmic and viscous heating in the accretion disc is extracted by 
a hot jet that is compatible with the result in CK04. In order 
to illustrate the thermal effect on the outflow, we have plotted 
the temperature vertical profile in Fig.|6]along a radius located 
at 1 AU from the axis. In this plot, the temperature increases in 
the disc corona before reaching its maximum T = 10 5 K and 
remaining constant. 

To study the time evolution of both accretion and ejection phe- 
nomena in the accretion disc and around the star, we analyze 
the accretion and ejection mass loss rate in both components. 
As in CK04 we draw the time evolution in Fig. of the mass 
loss rate Mj etl D in the disc-driven jet normalized to the accreted 
mass rate M A j at the inner radius Rj = I. We also display M A j 
normalized to the fixed mass accreted M a .e at the external ra- 
dius of our accretion disc at Rj - 40. Similar to CK04, we 
observe a strong increase, on the accretion rate in the inner part 
with time (Fig.0. This behavior is related to the extraction of 
the rotational energy of the accretion disc by the magnetic field. 
Indeed the creation of the toroidal component of the magnetic 
field in the disc brakes the disc matter so that the centrifugal 
force decreases, leading to an enhanced accretion motion. The 
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Fig. 8. Plot of the various forces Left: along a given magnetic field line anchored in the accretion disc, Right: along a streamline 
anchored to the stellar corona. These plots show the various forces accelerating the flow: is the centrifugal force, the 
magnetic one, fp the pressure gradient, and /g the gravitational force. 



mass flux associated with the disc-driven jet slowly increases 
to reach 18% of the accreted mass rate at the inner radius and 
contributes to 98% of the total mass-loss rate of the outflow. In 
fact, in this simulation the mass-loss rate from the central object 
is constant (1CT 9 M Q /yr), while the inner accretion rate reaches 
10~ 6 M Q /yr and the disc-drivenjet mass rate 10~ 7 M G /yr. Hence 
the stellar outflow does not affect the overall structure of the 
outflow much. This is confirmed by the shape of the outflow 
since it reaches a very similar aspect to the one obtain in CK04 
or in the previous simulation without a stellar jet, i.e. a jet con- 
fined within 20 inner disc radius. 

In order to analyze this accretion-ejection engine, we have to 
identify the forces responsible for the establishment of a steady 
accretion motion in equilibrium with a continuous emission of 
matter at the surface of the accretion disc. Furthermore we have 
to look at the collimation of the outflow and its interaction with 
the stellar wind. We show in Fig.[8]the various forces parallel 
to a given magnetic field line anchored in the disc and a flow 
streamline anchored to the central object. The various forces 
working along and across the field lines are 
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Fig. 9. Plot of the transverse forces as a function of R at a given 
altitude Z = 75. It shows, across a given cross section of the jet, 
the collimation processes acting in the stellar and disc com- 
ponents of the jet for the simulation with a stellar mass loss 
M s = lO _9 M /yr. is the centrifugal force, fp the pressure 
gradient and f£ the gravitational force, the Lorentz force 
due to the toroidal component of the magnetic field, while 
corresponds to the poloidal component. 



Ju = Ju,+fu e = (JxB)-e 



(39) 



where / p 5 '", /*'*, J**, and correspond respectively 
to the pressure gradient, the gravitational, the centrifugal, and 
the magnetic forces, one induced by the poloidal magnetic field 
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Fig. 10. The vertical variations of specific energies, normalized 
to the maximum kinetic energy flux, along a streamline in the 
stellar wind with a rate M = 10~ 9 M Q jyr. This plot clearly illus- 
trates the thermal acceleration provided by the enthalpy to the 
stellar wind material. The enthalpy gradient is here sustained 
by local turbulent Ohmic heating representing, in our simula- 
tion, 35% of the energy released by accretion (see Sect. 13.2.21 



component and the other by the toroidal one. The indexes *, " 
denote forces parallel and perpendicular to the poloidal mag- 
netic field line, respectively. We have defined the unit vector 
e s that corresponds to B p /\B p \ in the left panel of Fig. [8] and 
to Vp/|V ; ,| in the right panel of Fig. [8] while the perpendicular 
unit vector is such that e n ■ e s = in Fig.|5] As the stellar wind 
is resistive, the flow streamlines do not have to be parallel to 
the poloidal magnetic field to reach a quasi steady-state, since 

VpXB p = l] m Jg + 0. 

3.2.1 . MHD disc-driven jets 

The disc-driven jet behaves as in CK04. In particular, we find 
that the mass acceleration encounters two different regimes. In 
Fig.[8]we see that the vertical outflow is lifted from the accre- 
tion disc by both the magneto-centrifugal force and the pres- 
sure gradient up to the Alfven surface. Beyond, the poloidal 
acceleration is mainly sustained by the pressure gradient asso- 
ciated with the toroidal component of the magnetic field. Inside 
the resistive accretion disc, the toroidal component of the mag- 
netic field increases because of the differential rotation of the 
disc v m dB g /ds ~ J. dsrB p ■ VQ. Conversely, outside the disc, 
the advection of the toroidal magnetic field balanc es the dif- 
ferential rotation V ± ( Bg V „ ) = V(B p Q) exactly l lFerreirall 19971 
ICasse & Ferreiral2000al) . This change induces a decrease in Bg 
outside the accretion disc. As shown in CK04, this configura- 
tion allows matter below the disc surface to be pinched and to 
remain in an accretion regime, while beyond the disc surface, 
the change of sign of enables acceleration of mass along 



the magnetic field lines (cf Fig. [S}. This change in the mag- 
netic poloidal force also corresponds to a change in sign of the 
magnetic torque ((/ x B) ■ B p = -(J x B) ■ Bg) leading to the 
transformation of the MHD Poynting flux generated by the disc 
into kinetic energy of the jet material. 

The cylindrical collimation of the external outflow is induced 
by the pressure gradient of the poloidal component of the mag- 
netic field (f^ in Fig. [9}. In fact, the magnetic field in the disc- 
driven jet undergoes an expansion that induces a decrease in 
the poloidal magnetic field in the jet compared to the outer re- 
gion (Fig. Conversely the pressure gradient of the toroidal 
magnetic field (f^ ) acts to decollimate the outflow because the 
value of Bg is low outside the outflow, and its absolute value de- 
creases between the massive part of the disc-driven outflow and 
the outer medium. The magnetic field lines in the massive part 
of the outflow are anchored to the inner part of the accretion 
disc and extract more angular momentum than the magnetic 
lines in the outer medium. The inner part of the jet is, on the 
other hand, collimated by the toroidal pinching force. 

3.2.2. Stellar wind embedded in a disc-driven jet 

The stellar wind undergoes a thermal acceleration, as long as 
the shape of the jet does not become cylindrical. Incontrast 
to the disc-driven jet, the magneto-centrifugal force remains 
weak along the stellar wind flow. The opening angle between 
magnetic field lines that emerge from the sink region remains 
weak. In fact the stellar outflow starts to be collimated by the 
surrounding hollow jet induced by the accretio n disc. This ex- 
plains the difference from a model of iMatt & Balickl d20Q4 
where the stellar wind is the only outflow and which is prone 
to a dipolar expansion. In our simulation, we self-consistently 
describe the acceleration of the inner jet in addition to its colli- 
matio n, something that can b e compared with analytical mod- 
eling (Sa utv et a ■2002112004 . Part of both the thermal energy 
deposited at the surface of the corona near the polar axis and 
the energy deposited by the turbulence in the stellar wind is 
transformed into kinetic energy (Fig. 1101 along the streamline. 
We have estimated the amount of Ohmic heating released in the 
stellar outflow in the context of our simulation, namely, 



Ohmic 



(T] m J 2 - B ■ (V X TJ m J))dV 



(40) 



where V stands for the stellar outflow volume. The volume- 
integrated heating represents 35% of the energy released by 
accretion. However, the streamlines in the stellar wind are sub- 
ject to a larger expansion (relatives to the jet cylindrical radius 
at the Alfven surface) than the streamlines anchored to the ac- 
cretion disc. Therefore, the angular momentum conservation 
induces a decrease in Vg and Bg. In the asymptotic region, the 
magnetic field lines become almost radial so that the projec- 
tion of the magnetic pressure gradients along the magnetic field 
lines become positive. The flow undergoes magnetic and ther- 
mal acceleration in this region. 

The collimation of the inner part of the jet is induced by the 
thermal pressure plus the pinching of the toroidal component 
of the magnetic field. They balance the centrifugal force and 
the pressure of the poloidal component of the magnetic field 
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Fig. 11. Same plots as in Fig.[S]but for a simulation where the stellar mass ejection rate is M — 10 1 M Q /yr. 



(Fig- E}- Besides that, the simulation shows that the inner por- 
tion of the stellar jet has a deep in density (Fig@} and a peak of 
velocity around the axis ("Fig. ll3t . Thus, all these facts suggest 
that the very inner part of the outflow, the so-calle d spine jet, 
behav e as a meridionally-self-similar jet as in lSautv et all f2002. 
120041) . We observe that this is a kind of "hollow" stellar jet 
"thermally" driven and both magnetically and thermally con- 
fined inside the "hollow" disc jet. This result is to be compared 
to the analysis of CK04 where it is shown that the external disc 
jet is, partially at least, behaving as a radially self-similar one. 
In Fig.^] we display temperature isocontours within a small 
area around the sink region. Thanks to this plot, we can see 
that the magnetic lever arm associated with the various outflow 
components is different. In fact, the disc-driven jet exhibits a 
magnetic lever arm (related to the ratio of the Alfven radius 
to the magnetic-field line foot-point radius) varying approxi- 
mately from 9 to 25, while the magnetic lever arm associated 
with the stellar wind ranges from near the axis to several tens, 
if one considers the foot-point of magnetospheric magnetic- 
field line to be anchored to the star. This last magnetic lever 
arm value may not be very reliable since we have imposed the 
size of the sink and thus influenced the radial extension of the 
magnetospheric outflow near the sink. 



3.2.3. Ideal MHD two-components outflows 

An alternative to the presence of turbulence in the stellar wind 
would be to have a mechanism acting near the star corona and 
transforming a part of the accretion energy into thermal energy 
(see for instance Matt & Pudritz 2005 and reference therein). 
We performed several simulations without turbulence inside 
the stellar wind (a w = 0) and with the value of 5 E regularly 
increased (and thus the amount of thermal energy at the top of 
the sink). We find that, in our simulations, the maximal value 
of 6 S is around 10~ 3 . Beyond that value, the pressure above the 
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Fig. 12. Same density plot as in Fig.|3]but for a smaller zone 
around the sink. The three critical surfaces are represented 
by dark lines (slow-magnetosonic), dashed lines (Alfven), and 
white lines (fast-magnetsonic). The size of the sink region is 
Rj = 0.1 AU and the stellar mass is IM Q . The sink region is 
represented by a black square in the center 



sink is so high that it disrupts the accretion disc structure and 
prevents the launching of the disc-driven jet. It is noteworthy 
that this value of S e is linked to very high value of thermal en- 
ergy released in the star corona. 

Indeed, the reader has to keep in mind that the top bound- 
ary of the sink is quite far from the stellar surface (typically 
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Fig. 13. The transverse variation of the vertical velocity at Z = 8AU, for the simulation with stellar mass loss M = lO- 9 M Q /yr, 
M, = M Q , and Ri = Q.1AU, Left: without any non-ideal effects in the stellar wind but with a large amount of thermal energy 
deposited at the base of the stellar outflow (S e = 10~ 3 ). Right: With turbulent viscosity and resistivity in the stellar wind and 
a small amount of thermal energy at the base (8 e = 10~ 5 ). The turbulence allows a better thermal stellar mass acceleration as 
velocity becomes higer near the polar axis. 
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20 stellar radii) so that if the flow undergoes a spherical expan- 
sion with a constant thermal energy flux, the thermal energy 
deposited in the corona would represent S e {Ri/R*) 2 of the en- 
ergy released by accretion. This amount of thermal energy may 
then represent a significant fraction of the accretion energy. In 
our calculations done without stellar outflow turbulence, we 
noticed that the structures fulfilling observational constraints 
coincide with the highest value of 5 S allowing disc-driven jet 
launching (typically 10" 3 ). 

The resulting two-component outflow is very similar to 
simulations done with turbulence in the stellar wind (and a 
very low 6 S ), except the terminal velocity of the stellar out- 
flow as shown in Fig.^]where we display the vertical velocity 
of matter along a radial direction located at Z = 807?,. In this 
figure we can clearly see that the stellar wind prone to turbulent 
heating is faster than the ideal MHD stellar wind. The poloidal 
mass acceleration in this zone is very sensitive to thermal heat- 
ing since magneto-centrifugal is vanishing here. A continuous 
heating, as generated by Ohmic heating, seems more efficient 
for accelerating the mass since it is "refilling" the thermal en- 
ergy reservoir available for acceleration along the flow. 



3.3. Massive stellar winds vs. sun-like mass-loss rate 
wind effects on a disc-driven jet 

In the simulations presented so far, we have seen that winds 
with a mass-loss rate similar to the Sun (up to 10~ 9 M Q /yr) do 
not greatly influence the disc outflow since their general be- 
havior remains similar. However in the case of a massive stel- 
lar jet, the inner wind may strongly influence the outflow as 
can be seen in a new simulation performed for a stellar-wind 
mass-loss rate set to 10~ 7 M Q /yr (Fig. I14> . The radial stellar 
wind strongly compresses the magnetic field anchored in the 
accretion disc. The enhanced magnetic field bending (even in 
the external part of the accretion disc R > 30) leads to an in- 
crease in the magnetic pinching in an extended region of the 
disc 1 < R < 30. Thus the outflow is launched from this whole 



region since the Blandford & Payne criterion is fulfilled ev- 
erywhere (Blandford & Pavne 1982). Indeed the magnetic field 
becomes dynamically dominant in the disc corona of this re- 
gion. The magnetic bending larger than 30° from the vertical 
direction leads to a centrifugal force and a thermal gradient 
pressure that is more efficient for launching the outflow from 
the disc (Fig. II 1> as can be seen in the jet mass loss that reached 
0.5 of the accretion rate in the inner part (Fig.ll5>. However, the 
amplitude of the magnetic force e p ■ (J x Bo) projected along 
the poloidal direction becomes weak since the projection along 
e p of the magnetic pinching force increases with the expansion 
of the magnetic field (see Fi g . 1 1 1 1 . 

The angular momentum carried away by the stellar outflow 
now represents 5% of the accreted angular momentum at the 
inner radius of the accretion disc. Regarding the acceleration 
of the outflow, we can distinguish two regions: an internal 
one corresponding to the contribution from the stellar out- 
flow and an external one coming from the disc-driven jet. This 
last component reaches velocities up to v z = 15 (Fig. 1 1 61 . 
The acceleration of this component is thermally and magneto- 
centrifugally driven, which is coherent with the larger radial 
expansion of the stream lines (see Fig. II It . In the inner stel- 
lar wind, the flow undergoes a weak expansion, and its veloc- 
ity does not exceed v z = 6.8. The acceleration of this com- 
ponent is achieved mainly via the thermal pressure, which is 
expected since the mass density is much higher than in previ- 
ous simulations where M - IQ^M^/yr. Let us note that the 
fast-magnetosonic Mach number remains higher than one for 
the whole outflow (Fig.ll6>. 

4. Summary and concluding remarks 

In this paper we present numerical simulations of the interac- 
tion between an accretion-ejection structure launching a disc- 
driven jet and a stellar wind emitted either from the central ob- 
ject and/or from its magnetosphere, in particular for the case 
of YSOs. In our framework, the accretion disc is near equipar- 
tition between thermal pressure and magnetic pressure where 
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Fig. 15. Plot of the temporal evolution of the ejection mass-loss rate from the accretion disc in a) the ratio of the stellar mass-loss 
rate to the ejection mass-loss rate from the accretion disc as a function of time in b) for the simulation with a stellar mass-loss 
rate of M = 10~ 7 M Q /yr. 





B 




p 


B A 




* J \ 



















1 




V P 


8 












6 






4 
















2 










Mfasr 



40 60 
R (R, ) 



O 20 



40 60 80 

R <R, ) 



40 60 
R (R. ) 



Fig. 16. The transverse variation of different physical quantities as magnetic-field components, velocity components, and the 
fast-magnetosonic Mach number at Z = 100, for the simulation with stellar mass loss M = 10~ 7 M G /yr 



turbulence is believed to occur. This turbulence is character- 
ized by a time and space-dependent anomalous resistivity and 
viscosity set by using an a description. The origin of the turbu- 
lence is still unknown but is not likely to arise from magneto- 
rotational instability since an equipartition disc is inconsis- 
tent w ith the development of such instability ( Qgilvi e~& Livid 
120011) . 

The properties of both the accretion disc and outflow were in- 
vestigated in this paper. In a first stage, we analyzed the con- 
tribution of the various hydrodynamical and magnetohydrody- 
namical mechanisms to the angular momentum transport in the 
thin accretion disc including, for the first time, both anomalous 
viscosity and resistivity, with a magnetic Prandtl number equal 
to unity. We demonstrated that the MHD Poynting flux associ- 
ated with the disc-driven jet plays a major role in the removal 
of the angular momentum from the thin accretion disc. The an- 
gular momentum radial transport provided by the anomalous 
viscosity inside the disc remains weak and contributes only 
1% of the total angular momentum transport (this value agrees 
with analytical estimates depending on disc thickness and a 
value). This is consistent with the viscous torque depending 
upon the radial derivative of angular velocity, while the mag- 
netic torque is mainly controlled by the vertical derivative of 



the toroidal component of the magnetic field. In a thin accretion 
disc, the toroidal magnetic field varies from zero to an equipar- 
tition value on a disc scale height eR, e <K 1, leading to a much 
more efficient extraction of the rotational energy from the mag- 
netic torque into the MHD Poynting energy flux feeding the jet. 
Basically, with our simulation and despite the disc viscosity, we 
have reproduced very similar results to those obtained by CK04 
where the viscous torque was neglected. 
In the second stage of the present paper, we studied the effects 
induced by the launching of a stellar wind inside the hollow jet 
arising from the accretion disc. The stellar outflow is thermally 
driven by the turbulent viscous and resistive stresses in a region 
close to the polar axis. In order to mimic the non-ideal effects 
believed to occur inside the stellar wind, we prescribed anoma- 
lous viscosity and resistivity in the wind region, leading to a 
turbulent heating of the plasma near the polar axis. 
We performed various simulations using different stellar mass 
ejection rates from the central objects. These stellar ejection- 
mass rates range from a Sun-like star (M = 10~ 9 M G /yr) to O 
and B type stars (M = 10~ 7 M G /yr). The influence of the stellar 
wind on the dynamics and the structure of jet and the accretion- 
ejection structure around the stellar object gets stronger with 
higher stellar mass ejection rates. As an example in the simu- 
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Fig. 14. Same figure as in Fig. [2 but with a non-ideal stellar 
wind emitted from the inner region with an ejection mass rate 
M - 10~ 7 M Q /yr. The outflow structure is substantially mod- 
ified by the presence of the stellar outflow since its radial ex- 
tension is two times larger than in the case with no, or weak, 
stellar outflow. The size of the sink region is = 0.1 AU and 
the stellar mass 1M . 



lations where stellar mass ejection rates are M = 10~ 9 M Q /yr, 
we obtained a very similar disc-driven jet to the one in CK04. 
The only difference lies in the presence of the internal, fast, hot 
plasma coming from the central object. In this simulation, the 
ejection-mass rate in the disc-driven jet is similar to the one in 
CK04, on the order of 15% of the inner disc accretion mass 
rate, while the stellar ejection mass rate represents 1 % of the 
total mass loss in the outflow. 

In our simulations, the collimation of the stellar outflow 
takes place once the jet from the accretion disc is launched 
and has reached a significant spatial extension. Its collima- 
tion is provided mainly by pinching the toroidal magnetic field 
in equilibrium with the thermal pressure gradient. Conversely 
the collimation of the jet from the accretion disc is induced 
by the poloidal magnetic field pressure gradient balancing the 
centrifugal force. Furthermore, in all the simulation the stellar 
wind keeps having a more or less conical expansion up to the 
asymptotic region where the disc-driven jet acts to collimate 
the stellar flow into a cylindrical flow. 



These important results are sel f-consistently obtained, 
in contrast to simulations of Bogov alov & Tsinganosl |2005) 
where a relativistic wind was collimated by a jet but which 
did not consider either the accretion disc, the jet launching, 
or the stellar wind acceleration. Indeed in our model, we self- 
consistent describe the launching and the collimation of disc 
and stellar wind. In particular, our simulations completely de- 
scribe both the stellar flow acceleration (the stellar flow is in- 
jected with sub-Alfvenic velocity) and the launching mecha- 
nism of the jet from the accretion disc. 

We have shown that the contribution of non-ideal MHD 
mechanisms in the acceleration of the stellar outflow can be 
significant since turning on this dissipative mechanism leads, 
for instance, to higher terminal velocities of the stellar jet- 
collimated flow. Our prescription of these dissipative mecha- 
nisms is of course subject to improvements, but our goal was 
to show that they enable an increase in the efficiency of both 
the thermal and magnetic acceleration of the stellar wind. The 
turbulence may be produced by the interaction of the stellar 
wind with the disc-driven jet. Moreover, as in the solar wind, 
the turbulence in the stellar wind may also have a stellar origin 
and/or a possible connection to the accretion occurring near the 
surface of the star. In this scenario, a part of the energy released 
by accretion is carried away in the wind by outwardly propagat- 
ing Alfven waves inducing turbulence. This scenario is a vari- 
ant of models where a significant part of the accretion energy 
is converted into thermal energy in the star corona (see Matt 
& Pudritz 2005 and references therein). By performing simu- 
lations with no stellar wind turbulence but with a large amount 
of thermal energy at the base of the wind, we found quite sim- 
ilar result except for the velocity field, the resistive continuous 
heating of the stellar wind being more efficient in providing 
higher velocity. It is noteworthy that, similarl to models de- 
positing thermal energy near the stellar corona, the amount of 
energy released by the turbulent heating is a significant fraction 
of the accretion energy (in the particular case of our simula- 
tions, it represents near 35% of the accretion energy). Note also 
that similar double-layer jets were found by Koide and collab- 
orators in various simulations (e.g. Koide et al. 1998; Koide et 
al. 1999) but those simulations were devoted to rapidly variable 
jets (with only a few disc rotations) and not to steady structures. 
We performed simulations with higher stellar ejection-mass 
rate, typically with M = 10~ 7 M G /yr (compatible with O-B 
type stars). The increase in the stellar mass-loss rate induces 
a faster and larger expansion of the jet. Indeed the enhanced 
pressure provoked by the stellar wind tends to bend the disc 
magnetic field lines over a larger radial extension, leading to a 
larger disc-driven jet. The corresponding disc-driven jet mass- 
ejection rate is much higher than in previous simulations since 
it reaches 50% of the disc inner accretion rate (stellar ejection- 
mass rate on the order of 10%). The simulations give a quan- 
titative threshold beyond which the stellar jet gives a signifi- 
cant extra expansion of the disc jet. Typically a mass loss rate 
from the star on the order of M — 10~ 7 M e /yr gives a factor 
two in the radial expansion of the disc-driven jet. Although the 
total jet remains small in cross section, as in CK04, and com- 
pared to self-similar disc wind models. Moreover the stellar jet 
play a very important role, also, in the formation of the disc 
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wind. We verify that, part of the external disc wind may look 
quasi radially-self-similar in nature, the most inner part of the 
stellar wind is quasi meridionally self-similar. In our work we 
neglected all radiative losses coming from the central star or 
from the plasma itself. The implementation of these terms and 
the study of their impact on the outflow structure is postponed 
to future work. 
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